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Abstract. Next generation neutrino oscillation experiments will need a quantitative understanding of neutrino-nucleus in- 
teraction far better than ever. Kinematics covered by the relevant neutrino-nucleus interaction spans wide region, from the 
quasi-elastic, through the resonance region, to the deeply inelastic scattering region. The neutrino-nucleus interaction in each 
region has quite different characteristics. Obviously, it is essential to combine different expertise to construct a unified model 
that covers all the kinematical region of the neutrino-nucleus interaction. Recently, several experimentalists and theorists got 
together to form a collaboration to tackle this problem. In this contribution, we report the collaboration's recent activity and a 
goal in near future. 
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INTRODUCTION 

The breakthrough measurement of non-zero 613 is changing the neutrino oscillation experiments into a new phase. 
The next target will be the leptonic CP violation, and the mass hierarchy of the neutrino. For aiming at this target, we 
need not only better controlled, higher statistics experiments but also a quantitative understanding of neutrino-nucleus 
interaction at the level of 5% or better The relevant neutii no-nucleus interaction covers a wide kinematical region, as 
described in Fig. 1 (left). From the low to high energy side, characteiistics of the neutrino-nucleus interaction changes: 
from the quasi-elastic (QE), through the resonance (RES) region, to the deeply inelastic scattering (DIS) region. 
Contribution from each reaction mechanism to the neutrino-nucleus interaction is displayed in Fig. 1 (right). There, 
CC (NC) stands for the charged-cuiTent (neutral-cuiTent). The single n processes are associated with the resonance 
excitation. A long baseline accelerator neutrino experiment like T2K utilizes the neutrino-nucleus interaction at the 
left-bottom corner of Fig. 1 (left) where QE and the single pion production due to the A resonance excitation are main 
mechanism. Meanwhile, an atmospheric neutrino experiment expects to find effects of CP-violation/mass hierarchy 
in the higher resonance and DIS regions, as indicated in the figure. Obviously, it is essential to combine different 
expertise to construct a unified model covering all the kinematical region for the neutrino-nucleus interaction. 

Recently, several experimentalists and theorists got together to form a collaboration at J-PARC branch of KEK 
Theory Center to tackle this problem [1]. Our strategy is to develop and combine baseline models for QE, RES and 
DIS. In this contribution, we repoit the cuiTent status of the baseline models for QE, RES and DIS, and discuss a future 
perspective. 



QUASI-ELASTIC REGION 



For the neutrino energy between 0.1 GeV and 1 GeV, the QE process dominates the neutrino-nucleus interaction, 
as shown in Fig. 1 (right). At T2K energy, this QE process is the most dominant and important. Since the neutrino 




FIGURE 1. (Left) Kinematical region of neutrino-nucleus interaction relevant to the next-generation neutrino oscillation exper- 
iments. The neutrino energy and squared four-momentum transfer are denoted by v and 2^, respectively. (Right) Contribution of 
each reaction mechanism to the neutrino-nucleus interaction. The symbol CT denotes the cross section. 



oscillations depend on the energy of the neutrino, the neutrino oscillation experiments must reconstruct the neutrino 
energy by selecting events from the CC QE processes. 

At the summary session of the Nulntl2 Workshop, we generally agreed that the spectral function of nuclei should 
be used to calculate QE process and pion production in the neutrino-nucleus interactions, rather than a simple uniform 
nucleon momentum distribution based on the Fermi-Gas model, in order to achieve better accuracy of the calculations. 
The nuclear spectral functions P{p,E) is the probability of removing a nucleon of momentum (p) from ground state 
of the nucleus leaving the residual nucleus with excitation energy (£). The nuclear spectral function naturally enables 
us not only to calculate both CC and NC QE cross section, but also to estimate a spectroscopic factor (Sp) that is 
the probability of the residual nucleus in a specific shell state. Basic calculations with the spectral function approach 
are fully described in Ref.[2]. Comparisons of electron-nucleus scattering data with the spectral function approach 
and with the Fermi-Gas model are also shown in Ref.[2]. This approach is also useful to calculate CC and NC QE 
processes followed by a 7 emission [3] for which the branching fraction is significant 40%). 

A remaining issue is to better understand effects of 2p-2h (2-particle-2-hole) configuration and/or MFC (meson- 
exchange current) in the QE process. This subject has been raised by several recent experiments on the QE cross 
section measurement [4]. The data seem to indicate a considerable contribution from 2p-2h and/or MFC effects, and 
have motivated active theoretical studies [5]. 



RESONANCE REGION 

The neutrino-nucleon interaction in the resonance region is a multi-channel reaction where not only single pion 
production but also two-pion production has a comparable contribution above the A(1232). Tj and kaon productions 
can also happen with a smaller probability. In order to deal with this kind of multi-channel reaction, an ideal approach 
is to develop a unitary coupled-channels model. In the following subsections, we discuss the unitary dynamical 
coupled-channels (DCC) model, and its extension to forward neutrino-induced meson productions. 

Dynamical coupled-channels model 

In the DCC model [6, 7, 8], we consider 8 channels: yN ,nN ,r\N ,7tA,pN ,aN ,KA.,KT.. We solve a coupled-channel 
Lippmann-Schwinger equation that contains meson-exchange potentials and bare A^* excitation mechanisms, thereby 
obtaining unitary reaction amplitudes. We analyzed n{y)N — > nN ,r]N ^KA^KL reactions data simultaneously up to 
W = 2.1 GeV {W : total energy). The analysis includes fitting about 20,000 data points. The DCC model gives a 



CC v-proton / CC v-neutron CC v-neutron / CC v-proton 




W (GeV) W (GeV) 



FIGURE 2. The structure function F^iQ^ = 0) for the neutrino-induced meson productions from the DCC model. The solid (red), 
dashed (purple), dash-dotted (green), two-dotted dash (blue), and two-dash dotted (orange) curves are for the TlN, KTtN, T\N, KA 
and fCL reactions, respectively. The sum of them is given by the thick solid (black) curve. 



reasonable description of meson production data in the resonance region, thus providing a good basis with which we 
proceed to the neutrino reactions. 

Forward neutrino-induced meson productions based on PCAC 

In the forward limit, the matrix element for a neutrino-induced meson production off the nucleon is given by the 
divergence of the axial current, accompanied by some kinematical factor With the Partially Conserved Axial Current 
(PCAC) hypothesis, we can relate the matrix element with those of the pion-induced meson production. Thus with 
TzN X {X = 7zN, nnN, rjNjKAjKT.) amplitude, we can calculate the forward neutrino-induced meson productions. 
Recently, We have done such an application with the DCC model for 7zN X [9]. We present the result in Fig. 2 
where the structure function F2{Q^ = 0) is shown as a function of the total energy W of the hadronic final state. 
While nN production is the dominant process up to W = 1 .5 GeV, above that energy, the TtTtN production becomes 
comparable to nN, showing the importance of the nnN channel in the resonance region above A(1232). Other meson 
productions, T7A^, KA, and KT. reactions have much smaller contribution. We remark that this is the first prediction of 
the neutrino-induced nnN, rjN, KY production rates based on a model that has been extensively tested by data. 

Future plan 

The cuiTent model for the neutrino reaction is limited to only the forward kinematics. In near future, we extend the 
model so that we can analyze the neutrino reaction of all kinematics in the resonance region. For that, we need to 
develop a dynamical model for both vector and axial currents. Regarding the vector cuiTent, we akeady have some 
part by having analyzed data for meson photoproductions on the proton. We will also analyze the neutron target data, 
so that we can decompose the vector current into isovector and isoscalar pieces; this decomposition is necessary for 
calculating the neutrino processes. The 2^-dependence of the vector current couplings will be fixed by analyzing pion 
electroproduction data. For the axial current, we develop a model consisting of non-resonant mechanisms and also bare 
A^* excitation mechanisms. The axial current couplings of A^-A^* transition can be determined from the corresponding 
nNN* coupling using the PCAC relation. The 2^-dependence of the couplings can be fixed using available data of 
the neutrino-induced single pion production. In this way, we obtain all ingredients to calculate the neutrino-induced 
meson productions for finite in the resonance region. 

DEEP INELASTIC SCATTERING REGION 

The deep inelastic region of neutrino-nucleon interactions corresponds to the kinematical conditions of > 1 GeV^ 
and > 4 GeV^, and the cross sections or structure functions are described by the parton model. The neutrino- 



nucleon cross section via charged current is expressed by the three structure functions Fi, F2, and Fy. 
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where ± indicates + (— ) for v (v), x is the Bjorken seating variable x = ^ / {2p ■ q) with the nucleon momentum 
p, the H^-boson momentum q, and = —q^, Mn and Mw are the nucleon and W-boson masses, y is defined by 
y = P ■ q/{p ■ k) where k is the neutrino momentum, and E is the neutrino-beam energy. In the leading twist and 
the leading order (LO) of the running coupUng constant a^, the structure functions are expressed in terms of parton 
distribution functions (PDFs): 



2xFi=F2, F2^ = 2x{d + s + a + c), F2'' = 2x{u + c + d + s), 
2x{d + s — u — c), xF^^ = 2x{u + c — d — s). 



xf;^ 



(2) 



The next-to-leading order (NLO) expressions are calculated by the convolution integral of the above PDFs with 
NLO coefficient functions. For describing neutrino-nucleus cross sections, the nucleonic structure functions fi,2,3 
should be replaced by the corresponding nuclear ones 3' which include nuclear corrections. The nuclear structure 
functions are expressed by nuclear PDFs (NPDFs) ff (x) for the parton type / in the same way as Eq. (2). They contain 
nuclear correction factors Wi{x,A,Z), which indicate deviations from a simple addition of proton (p) and neutron (n) 
contributions: 

ff{x,Ql) = Wi{x,A,Z) i [Zff{x,Ql) + {A-Z)ff{x,Ql)] , (3) 

where Z is the atomic number of a nucleus, A is the mass number, and Qq is the initial scale which is taken as = 1 
GeV2. The weight functions Wi{x,A,Z) are determined by a global analysis of experimental data on lepton-nucleus 
deep inelastic scattering (DIS) and Drell-Yan processes with nuclear targets [10]. 

Neutrino DIS data can be included in the global analysis for extracting the NPDFs. However, one should be careful 
that nuclear modification data such as F2/F2, where D indicates the deuteron, are not available in the neutrino 
scattering. Recently, there are some discussions on possible differences between nuclear modifications of charged- 
lepton and neutrino reactions. In the MINER vA project, structure functions of light nuclei will be measured, so that 
such an issue will be investigated independentiy. 

In Fig. 1 (left), we notice a region which does not belong to either the resonance region or the DIS one at 2 < 1 
GeV^ and > 4 GeV^. In order to provide a complete interaction model for various neutrino reactions, we need to 
describe such a region. At 0, there is a guiding principle like PCAC to connect a structure function to the pion 

scattering cross section. Together with the update of the NPDFs including new data, an appropriate description of this 
regions should be developed. These studies are currently in progress. 
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